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Detection of HCN Direct £-Type Transitions Probing Hot 
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ABSTRACT 

We report the detection of direct Gtype transitions towards the proto¬ 
planetary nebula CRL 618 during a study of vibrationally excited carbon chains. 
The J = 8,10,11,12,13,14 AJ = 0 transitions of HCN in its hrst excited bend¬ 
ing mode ^2 = 1 were detected in absorption against the continuum of the central 
H II region making use of the Effelsberg 100 m telescope and the Very Large Ar¬ 
ray. Additionally, the J = 9 direct Ctype transition was detected in emission 
presumably indicating a weak maser. All lines are blueshifted with respect to 
the systemic velocity of CRL 618 indicating that the lines originate from a hot, 
expanding circumstellar envelope. The HCN column density along the line of 
sight in front of the continuum is 2 x 10^®cm“^. 

Subject headings: line: identihcation — radiation mechanisms: non-thermal - 
radio lines: stars — stars: individual (CRL618) — stars: AGB and post-AGB 
— stars: mass loss and circumstellar matter 


1. Introduction 

Molecules in vibrationally excited states are unique tools to study hot molecular gas in 
interstellar and circumstellar environments. The vibrationally excited states can be pumped 
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by IR radiation in the dusty environment close to the exciting object (see e.g. Schilke et ah 
1992; Schilke et ah 2000 or Wyrowski et ah 1999 and references therein). 

A linear molecule with N atoms has 3N — 5 vibrational modes. For HCN, this yields 
four vibrational modes^, two of which are stretching modes (CN stretch ui, CH stretch z/ 3 ), 
and the H — C — N bending mode (z/ 2 ) which is doubly degenerate since the molecule is free 
to bend in two orthogonal planes. If the molecule is bending and rotating simultaneously, 
the degeneracy is lifted giving rise to a phenomenon denoted i-type doubling: The degenerate 
bending state can be regarded as having components of an additional angular momentum 
p = ih about the hgure axis with i = — 2,n — 4,..., —v. For a hrst excited bending 

mode i = +1, —1 causes the splitting of every rotational level into two sublevels. Now two 
different types of transitions with either AJ = ±1 or AJ = 0 may occur, the latter being 
denoted direct i-type transitions (cf. Figure 1). Since the splitting of the sublevels of a given 
J is rather small, the corresponding transitions occur at low frequencies. For a molecule in 
a hrst excited bending mode those are to hrst order given hj u = qJ{J 1), where q is the 
Atype doubling constant {q{ucN,v 2 =i) ~ 224MHz). 

In the laboratory, direct Atype transitions were hrst measured by Shulman and Townes 
(1950) for the molecules HCN and OCS. In space, rotational transitions of vibrationally 
excited HCN (AJ = ±1, V 2 = 1,2; “vibrational satellites”) were hrst detected by Ziurys 
and Turner (1986) toward Orion-KL and IRC+10216. Meanwhile, HCN was also detected in 
much higher vibrationally excited states (e.g. V 2 = 4, Schilke et ah 2000; Schilke and Menten 
2002). A detection of the J = 19 direct Atype transition of HCN {v 2 = 1) in IRC+10216 
has previously been reported by Cernicharo et al. (1996), although no details are given. 

We started a detailed study of vibrationally excited carbon chain molecules toward the 
protoplanetary nebula CRT 618 at radio and millimeter wavelengths to examine their radial 
distribution by measuring their vibrational temperatures. In this letter, we report the hrst 
detection of seven consecutive direct Atype transitions of HCN (^2 = 1) at cm wavelengths, 
commencing with J = 8 through 14. The results on the higher cyanopolyynes HC 3 N and 
HC 5 N are reported separately (Wyrowski et al. 2002; Thorwirth 2001). 


"‘For the sake of clarity it should be recalled that i' is used to denote a certain vibrational mode or the 
frequency whereas v denotes the vibrational quantum number. 
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2. Observations 

2.1. Effelsberg 100 m observations 

The observations were performed during 8 observing sessions from 1999 March to 2002 
January employing the 100 m radiotelescope of the Max-Planck-Institut fiir Radioastronomie 
at Effelsberg, Germany. As frontends the 0.65 cm, 1cm, 1.3 cm, and 1.9 cm HEMT receivers 
with typical receiver temperatures from 30 to 70 K, respectively, were used. The FWHM 
beam widths and the frequencies of the HCN lines are given in Table 1. Pointing was checked 
every one to two hours on CRL 618 itself or other appropriate radio sources resulting in an 
average pointing accuracy of 4". The intensity scale was established using continuum drift 
scans on W3(OH), NGC 7027 and 3G147 and comparing with the flux densities given by Ott 
et ah (1994) . Pointing scans on GRL618 were used to determine its flux density relative to 
these calibration sources (Table 1). 


2.2. CSO observations 

Spectra of the J = 4 — 3 line of HGN in its V 2 bending mode were obtained in 1999 
December with the 10.4 m telescope of the Galtech Submillimeter Observatory (GSO)^ with 
a system temperature of 600 K. The line was observed in the upper sideband and several 
different observing frequencies were used to avoid blending from lines in the lower sideband. 
The spectrometer and observing procedure are described by Menten and Young (1995). The 
GSO beamsize at 356 GHz is 26" and we assumed a beam efficiency of 78% (taken from the 
GSO webpage). 


2.3. VLA observations 

GRL 618 was observed in the J = 13 direct Gtype transition (cf. Table 1) with the 
Very Large Array (VLA®) in its B configuration, leading to an angular resolution of 0.12". 
At the time of the observations, 15 antennas were equipped with 0.7cm receivers. A total 
bandwidth of 6.25 MHz was observed with 128 channels and the spectral resolution was 


^The CSO is operated by the California Institute of Technology under funding from the National Science 
Foundation, Grant No. AST-9980846. 

®The VLA is operated by the National Astronomy Observatory, a facility of the National Science Foun¬ 
dation operated under cooperative agreement by Associated Universities, Inc. 
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48.8 kHz. The total time on source was 3.5 h. Regular observations of 0555+398, 3C48 and 
3C84 were used for amplitude, flux and bandpass calibration, respectively. The phase was 
selfcalibrated on the strong continuum emission from CRL618. 

The remaining non Q-band antennas were used to observe the J = 4 direct f'-type 
transition at 4488.48 MHz. No line was detected at an RMS noise level of 4 mJy in 0.8 km/s 
wide channels. The total continuum flux at this frequency is 26 ± 3 mJy. 


3. Results and Discussion 

Figure 2 shows the direct f'-type lines observed toward CRT 618 and the line parameters 
from Gaussian fits to the spectra are given in Table 1. The J = 8,10,11,12,13,14 direct 
Gtype transitions appear as absorption lines towards the continuum of the H ll region with 
line velocities of approximately vl = —27kms“^. The systemic velocity Vgys of CRT 618 is 
—24.2kms“^ (Wyrowski et ah 2002), hence all of the observed lines are blueshifed relative 
to Vsys indicating that the lines originate from a hot, expanding circumstellar envelope. To 
determine the physical conditions of the absorbing gas, we used a spherical LTE model 
of an expanding envelope, developed to interpret our observations of vibrationally excited 
HC 3 N derived by (Wyrowski et al. 2002). To £t the HCN lines we use the temperature, 
density and velocity structure of the expanding envelope, which fits the HC 3 N lines, and 
only vary the HCN abundance. Using a temperature of 560 K (see discussion below) a HCN 
column density of 2 x 10^® cm“^ is needed to cause the observed absorption. The resulting 
fits to the spectra are shown in Figure 3 together with a spectrum of the HC 3 N = 1 
J = 12 — 11 transition which has a similar upper energy (~ 1300 K) Since no continuum 
flux measurements were performed at 350 GHz, the HCN J = 4 — 3 spectra are shown in 
brightness temperature units, whereas for the other spectra the ratio of line to continuum 
temperature, which reduces calibration uncertainties, is shown. The deviation of the HCN 
J = 4 — 3 model from the observed spectrum could be due to pointing and/or focus errors: 
a pointing error of 10 " alone would explain the difference between observation and model 
and cannot be excluded. The model consists of power laws for temperature, density and 
velocity starting at an inner radius of 0 . 11 " and an Hll region within that radius, which fits 
the continuum measurements of CRL 618. The temperature at the inner radius is 560 K. 
In the model, the HCN J = 4 — 3 line is highly optically thick and mostly sensitive to the 
model temperature and the emitting size. The HCN direct Ctype lines, on the other hand, 
are optically thin and probe a combination of temperature and column density of the model. 
The best fit model has a HCN/HC 3 N abundance ratio of 3 to 6 dependent on the assumed 
population of the vibrational levels of HC 3 N which is consistent with the result obtained by 
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mid-infrared absorption measurements by Cernicharo et al. (2001). 

Figure 4 shows the results of the VLA observations. The total continuum flux density 
at 40 GHz, estimated from the flux on the shortest baselines, is 0.75 Jy with an uncertainty 
of 10%. The size of the continuum emission is 0.34 x 0.16", estimated from a Gaussian £t to 
the UV data. To increase the spectral sensitivity, every four channels were averaged together 
and a taper in the UV plane was applied, reducing the angular resolution to 0.34 x 0.31". The 
insert in Figure 4 shows the spectrum integrated over the indicated area. Line parameters 
of a Gaussian fit to the spectrum are given in Table 1. To image the HGN absorption the 
continuum was subtracted from the data using the channel ranges marked in the insert of 
Figure 4. The contours in Figure 4 show the HGN absorption averaged over the line. 

The observed 40 GHz continuum emission compares well with the results of Martin- 
Pintado et al. (1993, 1995) at 23 GHz. The HGN absorption falls into the same velocity 
range as the hot core (HG) component seen in ammonia by Martin-Pintado & Bachiller 
(1992). The HGN absorption is slightly shifted to the west from the center of the continuum, 
which was also observed for the HG ammonia component and the hot dense disk observed 
by Mart in-Pint ado et al. (1995). However, no absorption is observed at the ammonia broad 
absorption velocity of —50kms“^ which Martin-Pintado et al. (1993, 1995) interpreted as 
occurring from post-shocked clumps. Accordingly, the HGN absorption most likely originates 
from the same volume of gas as the ammonia hot core component and the dense disk. 


A weak J = 9 maser? 

In contrast to the J = 8,10,11,12,13,14 transitions appearing as absorption lines to¬ 
ward the continuum, the J = 9 transition is found in emission. Blending with an unknown 
line cannot be ruled out entirely, but to the best of our knowledge there is no known tran¬ 
sition of a different molecule, and no trace of a J = 9 absorption (which would modify the 
emission-profile of the blending line) is seen. Moreover, the (blueshifted) velocity and the 
line width correspond well to the velocities and linewidths of the J = 8,10,11,12,13,14 
transitions. In particular the linewidth argument is compelling, since other emission lines 
are much broader (e.g. the linewidth of the HG 7 N J = 21 — 20 transition is 27 km s“^, 
Martin-Pintado & Bachiller 1992). These facts suggest that the J = 9 transition is a weak 
maser amplifying the continuum. Its optical depth has a different sign, but is similar in value 
to the absorption lines. 

How can this maser be understood? A quantitative analysis would require quite detailed 
modeling of the pumping mechanism, which is not feasible because many of the transition 
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rates (in particular collision rates between vibrational states and within the vibrationally 
excited state) are not or only poorly known. However, one can argue as follows that this 
transition is easily perturbed or inverted: The splitting due to t'-type doubling is much 
smaller than the rotational splitting. Considering the rotational level system J = 8, J = 9 
and assuming that the excitation temperature of the J = 9e—8e is (and equal to 

while the excitation temperature for the / (upper) states is taken to be + AT 

(with AT <C Tg^x’®®), one can show that inversion occurs < 0) if following condition 

is met: 

AT > = 0.025 (1) 

^8f,8e 

In terms of occupation numbers, a difference AT = 2.5% means that for an excitation 
temperature of 560 K, the occupation of the J = 9/ level needs to be elevated in population 
by only 0.17% with respect to a completely thermalized distribution to invert the direct 
f'-type transition. If one considers the populations relative to the vibrational ground state, 
similar arguments can be made. A similar system is the ammonia molecule, where the 
inversion splitting is much smaller than the rotational splitting, and indeed, some masers 
are found for ammonia as well (see Schilke, Walmsley & Mauersberger 1991 or Madden et 
ah 1986). It has to be emphasized that the actual pumping mechanism is probably quite 
different for the ammonia masers, the similarity is that it is not difficult to produce masers 
in the inversion lines. 

Having argued that a small perturbation is sufficient, we still have to identify a possible 
cause for perturbations in our f'-type system. One good candidate is a line overlap. Such 
a mechanism has been invoked successfully for pumping of OH masers (e.g. Cesaroni & 
Walmsley 1991), although the details of the pumping mechanism are different. One possible 
candidate responsible for a perturbation in the J = 8,9 V 2 = 1 levels is an overlap of 
rotational lines of the vibrationally excited molecule: For the J = 8 — 7 and J = 9 — 8 
transitions of the V 2 = 1 and V 2 = 2 states the i components If and 2e are separated by 
only 11 and 58MHz^, respectively, corresponding to Av = 5kms“^ and Av = 22kms“^. 
Accordingly, the hrst overlap could occur locally whereas the second one could connect 
different parts of the envelope. A corresponding maser effect has been observed for SiO 
masers by Cernicharo et ah (1991, 1993). 

V = 8 - 7, W 2 = 2, ^ = 2e at 712361MHz; J = 8 - 7, uz = 1, f = 1/ at 712372MHz; J = 9 - 8, V 2 = 2, 
e = 2e at 801421MHz; J = 9 - 8, V 2 = 1, f = 1/ at 801363MHz 
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HCN Direct i-type transitions in High-Mass Star-Forming Regions 

Initiated by the results presented here, additional searches for direct t'-type transitions 
of HCN have been performed towards high-mass star forming regions. So far, we were able 
to detect the J = 9 transition toward Orion-KL and SgrB2(N) and the J = 9,10 transitions 
toward GIO.47-1-0.03 using the Effelsberg 100 m telescope (Thorwirth et ah, in preparation). 
Direct f'-type transitions are optically thin and hence can be used to derive reliable column 
density estimates of the hot gas component. Moreover, the VLA observations of CRL618 
presented here demonstrate that this special kind of transitions can be used to observe hot 
gas at high angular resolution and low frequencies. 


4. Summary 

Using the Effelsberg 100m telescope and the VLA, we have detected seven direct Atype 
transitions of HCN in its V 2 = 1 state. All lines appear in absorption against the embedded 
H II region, except for the J = 9 transition, which shows weak maser action. The observed 
line velocities agree well with those observed in vibrationally excited HC3N by Wyrowski et 
ah (2002), who used their extensive data to model the physical parameters of the hot, dense 
emission region. VLA observations of the J = 13 line reveal the emission region is compact 
and covers only the western part of the embedded H ll region, similarly to ammonia. Since in 
addition to CRL 618 we detected Atype HCN lines toward hot molecular cores in regions of 
high-mass star formation, they also represent an interesting new tool to study the immediate 
vicinity of young (proto-)stellar objects. 

Triggered by the radioastronomical observations presented here, a laboratory investi¬ 
gation of direct Atype transitions of HCN (^2 = 1) has been carried out in the Cologne 
laboratory covering rotational quantum numbers up to J = 35 at 278.7 GHz (Thorwirth et 
ah, in preparation) The complete analysis will be presented in a following paper. 


The present study was supported by the Deutsche Eorschungsgemeinschaft (DEG) via 
Grant SFB 494 and by special funding from the Ministry of Science of the Land Nordrhein- 
Westfalen. FW is supported by the National Science Foundation under Grant NSF AST- 
9981289. We would also like to thank an anonymous referee for valuable suggestions on the 
maser interpretation. 
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Table 1. Direct t'-type transitions of HCN (MHz) in its ^2 = 1 vibrational state observed 
in the present study, as well as corresponding upper energies (K), beam sizes 9mb {"), line 
to continuum ratios, line intensities II (mJy), line velocities vl (kms“^), line widths An 
(kms“^) and continuum flux densities (mJy). 


J 

Frequency ^ 

Eu 

^vnb 

-Tl/Tc 

h 

Vl 

An 

gcont 

4 

4488.48 

1067 


b 




26 (3) 

8 

16148.55 

1178 

52 

c 


-26.6(3) 

3.8(5) 


9 

20181.40 

1217 

42 


22 (4) 

-27.0 (3) 

4.7(7) 

425(60) 

10 

24660.31 

1259 

34 

0.11(3) 


-26.6(3) 

4.4(6) 

490 (90) 

11 

29584.66 

1306 

28 

0.05(1) 


-27.0 (2) 

4.8(6) 

d 

12 

34953.76 

1358 

24 

0.07(2) 


-27.4 (2) 

5.1(4) 

680(140) 

13 

40766.90 

1413 

0.33 

0.11(4) 


-27.3 (7) 

5.4(12) 

750(75) 

14 

47023.20 

1473 

20 

0.11(4) 


-26.7(4) 

4.9(11) 

840(85) 


^Laboratory frequencies taken from Maki (1974). 

'^No line detected. 

'^Line detected in absorption without valid calibration. 
'^No absolute flux calibrator available. 
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Fig. 2.— The J = 8, 9,10,11,12,13,14 AJ = 0 (direct Atype) transitions detected toward 
CRL618. The J = 8, 9,10,11,12,14 transitions were detected using the lOOm telescope and 
the J = 13 transitions using the VLA, respectively. The dashed line marks the line velocity 
vl = —27kms“^ indicating a signihcant blueshift from Vgys = —24.2kms“^ of CRL618. 
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Fig. 3.— Comparison of observed spectra of vibrationally excited HCN/HC 3 N with results 
of an expanding envelope model. 
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Fig. 4.— Results of the Q-band VLA observations. The continuum emission is shown in 
greyscale with steps of 10, 30, 50, 70, and 90% of the peak continuum flux (140 mJy/beam). 
The HCN absorption is shown in contours of —16.5, —27.5, —38.5 mJy/beam, averaged from 
Visr —24.7 to —30.3 km/s. The beamsizes are indicated in the lower left and are 0.13 by 
0.11"and 0.34 by 0.31", respectively. The insert in the upper left shows a spectrum of the 
emission integrated over the region within the dashed ellipse. 




